Towards Classification of Phase Transitions in Reaction— Diffusion Models 
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Equilibrium phase transitions are associated with rearrangements of minima of a (Lagrangian) 
potential. Treatment of nonequilibrium systems requires doubling of degrees of freedom, which may 
be often interpreted as a transition from the "coordinate" - to the "phase" -space representation. 
As a result, one has to deal with the Hamiltonian formulation of the field theory instead of the 
Lagrangian one. We suggest a classification scheme of phase transitions in reaction-diffusion models 
based on the topology of the phase portraits of corresponding Hamiltonians. In models with an 
absorbing state such a topology is fully determined by intersecting curves of zero "energy." We 
identify four families of topologically distinct classes of phase portraits stable upon renormalization 
group transformations. 

PACS numbers: 64.60.Ak, 05.40.-a, 64.60.Cn, 82.20.-w 



I. INTRODUCTION 

The last decade witnessed the rapid growth of interest 
in reaction-diffusion modela 1 T 2 T 3 T 4 T 5 T 6 T 7 1 8 1 9 1 10 T 11 1 12 i 13 1 14 i 15 . 
Such models are employed for a description of phenomena 
ranging from kinetics of chemical reactions to the evolu- 
tion of biological populations. The subject of particular 
interest is the description of dynamical phase transitions 
in reaction-diffusion kinetics. An important example is 
absorbing phase transitions. Upon such a transition the 
system goes from an active ( "living" ) phase to an absorb- 
ing ("dead") state with no escape from it. 

Grassberger and JansseniSiii realized that many of 
the absorbing-state transitions belong to the same uni- 
versality class as the directed percolation (DP) model. 
Since then, the DP universality class was extensively 
s tudied both analytica lly and numerically (see Refs. 0, 
I3ll8ll9l20l2ll22l23l24| ). The DP universality class is 
extremely robust. In fact, exceptions to the DP tran- 
sitions are rare. However, if the microscopic dynam- 
ics possesses additional symmetries, the universality 
class of the transition may be different. For exam- 
ple, parity conservation (PC) is known to drive the 
transition to a new distinct universality cl&sai^^^. 
Recently other possible universality classes have been 
studied 1 ! 27 ! 28 . 29 ! 30 ! 31 ! 32 . 33 ! 34 ! 35 ! 36 : 37 ! 38 ! 39 ! 40 . Among them 
is a pair contact process with diffusion (PCPD), the 
critical behavior of which has not yet been described 

ana 1...^ :a U --8.41.42.43.44.45.46.47.48.49.50.51.52.53.54.55.56 

Despite extensive accumulated knowledge, it seems 
that a guiding principle, allowing one to distinguish be- 
tween various types of the transitions, is still missing. 
The purpose of this paper is to suggest a simple scheme, 
providing at least, an educated guess regarding the uni- 
versality class of the reaction-diffusion model at hand. 
The scheme is based on the topology of the phase por- 
traits of the system's Hamiltonian. Before elaborating 
on it, let us briefly recall the corresponding strategy for 
equilibrium systems. 

An equilibrium system may be characterized by an ac- 



tion (energy) 



S = j d d x [D{Vq) 2 + V(q)] 



(1) 



written in terms of the order parameter q(x) (for sim- 
plicity we restrict ourselves to the one-component order 
parameter). The potential function V(q) encompasses in- 
formation about possible phase transitions. Specifically, 
one monitors the behavior of the minima of V(q) as a 
function of the control parameter to infer the existence 
and type of transition. For example, a wide class of mod- 
els may be described by a potential of the form: 



V(q) = hq + mq 2 + u q 4 



(2) 



For m < the system exhibits a first-order transition 
when h changes its sign (for d > 1) and the two minima 
interchange. In the symmetric case, h = 0, the system 
may undergo a second-order transition when the param- 
eter m is swept through zero, so a single minimum is 
split in two. Below the critical dimension d c — 4 this 
second-order transition is characterized by non-mean- 
field critical exponents. To find the exponents one typi- 
cally employs the renormalization group (RG) technique. 
The RG technique treats the coefficients m and u of the 
potential J2J as functions of the running spatial scale. 
For e = d c — d > the potential scales towards the non- 
trivial fixed point potential, V*(q) ^ 0, with m* ~ e 2 and 
u* ~ e. Notice that the action essentially specifies the 
Lagrangian field theory. 

One may ask whether a similar strategy exists for non- 
equilibrium phase transitions in reaction-diffusion sys- 
tems. To answer this question one has to recall that a 
description of nonequilibrium systems requires doubling 
of the degrees of freedom. There are various manifes- 
tations of this statement depending on the specific con- 
text. In quantum kinetics it is known as the Keldysh 
technique 5 ^. It employs time evolution along a closed 
contour, so one has to keep two copies of each field: one 
for the forward and another for the backward evolution. 
In a classical context the Martin-Siggia-Roso 58 i 59 i 60 i 61 i 62 
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method requires an additional set of fields to resolve the 
functional (5-functions of Langevin equations. Most im- 
portantly for the present subject, in reaction-diffusion 
kinetics the Doi-PelitiSiSi operator technique deals with 
the creation and annihilation operators for each reagent. 
Thus it employs two variables (or one complex field) for 
every real physical degree of freedom. (For a discussion 
of the connections between these techniques see, e.g., 
Refs. |6l6^ .1 

An important observation is that in all these exam- 
ples the two sets of fields (being properly transformed) 
may be considered as canonically conjugated variables. 
As a result, instead of the equilibrium order parameter 
q{x), one has to deal with the canonical pair: q(x,t) and 
p(x, t). Correspondingly a reaction-diffusion system may 
be described by the Hamiltonian action 

S= JdtJd d x{pd t q + DV P Vq-H R (p,q)} , (3) 

where the Hamiltonian H R (p, q) is determined by the set 
of reactions specific for a given model (see below) . 

Comparing Eqs. © and QJ, one notices that the re- 
action Hamiltonian Hji(p^q) plays a role similar to the 
effective potential V(q) in equilibrium statistical mechan- 
ics. Thus it is plausible that H R (p,q) may encode in- 
formation about possible nonequilibrium transitions in 
a way analogous to what V(q) does. Specifically, one 
wants to know what is the Hamiltonian analog of the 
potential minima, given by d q V = 0, in the Lagrangian 
formulation. The answer is that it is the classical equa- 
tions of motion: dtq = d p H R and dtp = —d q H R . One 
is looking, therefore, for a geometric way to picture the 
Hamiltonian equations of motion. We argue below that 
the way to do it (at least for the one-component models) 
is to consider the phase space trajectories in the (p, q) 
plane. Indeed, the classical equations of motion conserve 
"energy." Thus the phase-space trajectories are given by 
the curves Hr(j>, q) = const. Moreover, for systems with 
absorbing states the only trajectories which may inter- 
sect correspond to zero energy. As a result, the set of 
curves 

H R (p,q) = (4) 

determines entirely the topology of the phase space. 

The main message of this paper is that the curves 
specified by Eq. (gj and the corresponding topology of 
the phase portrait classify possible phase transitions in 
reaction-diffusion models. It is the web of the zero "en- 
ergy" trajectories which plays the role of minima of the 
V(q) potential in the equilibrium statistical mechanics. 
A topological rearrangement of this web as a function 
of the control parameter signals the existence of a phase 
transition. The corresponding topology is in one-to-one 
correspondence with the universality classes. 

Below we show that the number of distinct generic 
phase-space topologies (for one-component systems) is 
rather limited, indicating that all possible universality 



classes may be exhausted. Some of the topologies corre- 
spond to first-order transitions, while others to contin- 
uous ones [much the way potential @ contained both]. 
The latter class may develop nontrivial critical exponents 
below a certain critical dimension d c . The way to find 
these exponents is to follow the RG flow of constants of 
the Hamiltonian Hji(p, q) upon elimination of the small- 
scale fluctuations. While the Hamiltonian itself may be 
complicated, it is only the topology of the phase space 
that matters, not the specific shape of the curves. Any 
given topology may be modeled by a simple polynomial 
of q and p, again much the way the simple polynomial 
(J2J) suffices to describe many equilibrium systems. Thus 
one must follow only changes of the topology of the phase 
portrait upon RG transformations. One should also ver- 
ify that a given topology is stable upon RG transforma- 
tions, i.e., that it cannot be reduced to a more generic one 
by decimation. The resulting fixed-point topologies and 
corresponding fixed-point Hamiltonians H R (q,p) provide 
information about the universality classes. 

We found that the DP universality class (represented 
by the simplest triangular structure on the phase plane) 
serves as parent for a family of descending classes. Each 
subsequent class in the family is characterized by a min- 
imal number fc of particles needed to initiate reactions. 
We denote it as the fc-particle contact process with dif- 
fusion (kCPD). Here 1CPD is DP, while 2CPD is the 
PCPD (for a review see Ref. § and references therein). 
According to scaling analysis, above the upper critical 
dimension d c — 4/fc the kCPD's are characterized by 
mean-field critical exponents, e.g., (3 — 1 and v±_ = fc/2. 
For d < d c and k — 2, 3 we found that RG flows to a 
strong-coupling fixed-point that cannot be accessed in 
the e expansion (see also Refs. 8.67]). We also discuss 
the possible nature of the strong-coupling fixed point for 
k = 2. 

Similarly the PC universality class generates a family 
of classes, characterized by a minimal number fc, of in- 
coming particles required for all reactions. We call them 
fc-particle parity conserving (kPC). Their upper critical 
dimension is d c = 2/k. In addition to kCPD and kPC we 
identify two more families of universality classes. They 
both originate from reversible reactions which may go 
in both directions with different rates. We call them fc- 
particle reversible (kR) and fc-particle reversible parity 
conserving (kRPC). In both cases fc stays for a minimal 
number of incoming particles. Their critical dimensions 
are 2/fc and 2/(fc + 1) correspondingly. 

These four families seem to exhaust all possible con- 
tinuous transitions reachable by the tuning of a single 
control parameter and capable of exhibiting a non-mean- 
field behavior. This means that any phase portrait, topo- 
logically different from that of the four families, is unsta- 
ble upon renormalization. In the large-scale limit it flows 
towards one of the stable topologies. The latter are pro- 
tected by certain symmetries of the action against defor- 
mations introduced by the RG. The fluctuation-induced 
renormalization may be not effective if the space dimen- 
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sionality is sufficiently high. As a result, other topologies 
may appear to be stable and lead to different universality 
classes. However, by virtue of the ineffectiveness of fluc- 
tuations, such universality classes are bound to exhibit 
trivial mean-field behavior. 

The paper is organized as follows: in Sec. ^ we in- 
troduce reaction Hamiltonians and their phase portraits. 
Section Ifffl is devoted to models exhibiting the first- 
order transitions and discusses the topological structure 
of their phase portraits. In Sec. IIVI models of the DP 
universality class and its derivatives, kCPD's, are con- 
sidered. We demonstrate that the triangular topology of 
the phase portrait is the typical feature of all universality 
classes of this kind. In Sec. we consider the parity con- 
serving model and its generalizations, kPC's. The rect- 
angular phase portrait topology of the reversible reaction 
models, kR's and kRPC's, is discussed in Sec. lVfl Finally 
some conclusions and outlook are drawn in Sec. IVIII 



II. REACTION HAMILTONIANS AND PHASE 
PORTRAITS 

The standard way to introduce the "quantum" reac- 
tion Hamiltonian is by employing the creation and anni- 
hilation operator technique of Doi and PelitiiSiiiiSiSiS. 
Here we choose to follow a different, though completely 
equivalent, strategy 6 -. Consider a generic reaction that 
transforms k particles into m equivalent ones with the 
probability A: 



kA mA, 



(5) 



The corresponding master equation for the temporal evo- 
lution of the probability V n (t) of a configuration with n 
particles has the form 



d t Pn(t) = A 



k — m 
k 



-\-k—m 



(*) 



V n (t) 



(6) 

The two terms on the right hand side (RHS) represent the 
probabilities of "in" and "out" processes correspondingly. 
The master equation © is to be supplemented with an 
initial distribution, e.g., V n (0) = e~ n °n,Q /n\, the Poisson 
distribution with the mean value no, or V n (0) — S n . no , 
the fixed initial particle number. 

Let us define now a generating function as 



G(p,t)=J2p n Vn(t). 



(J) 



n=0 



The parameter p will play the role of the canonical mo- 
mentum; so far it has been introduced pure formally. The 
value p = 1 plays a special role. First, the conservation 
of probability implies the normalization condition 



G(l,t) = 1. 



(8) 



Second, the moments of V n (t) may be expressed through 
derivatives of the generating function at p = 1, e.g., 



i{t)) =Y,rtV n {t) = d p G{p,t)\ p=1 



(9) 



Knowing the generating function, one may find a proba- 
bility of having (integer) n particles at time t as V n (t) = 

fi£G(p,t)|p=oM 

In terms of the generating function the master equa- 
tion © may be identically rewritten as 



d t G = H R {p, q)G, 



(10) 



where the non-Hermitian normally ordered operator H R 
stays for 



Hr(pA) = 77 (P 



p k )q k . 



(11) 



Here we have introduced the "coordinate" operator 
. d 



^ dp ' 



(12) 



obeying the canonical commutation relation [q,p] = 1. 
Because of the obvious analogy with the Schrodinger 
equation, we shall refer to the operator H R as the Hamil- 
ton operator in the p representation. From the normal- 
ization condition, Eq. J8J, it follows that 



H R (p,q) 



0. 



(13) 



Any Hamiltonian derived from a probability conserving 
master equation necessarily satisfies this property. 

One can easily generalize this construction for the case 
where many reactions with rates \km take place at the 
same time. To this end one can simply algebraically add 
the corresponding partial Hamiltonians to obtain the full 
reaction Hamiltonian. If there is no particle production 
from the vacuum, i.e., k ^ for any m-the empty state 
with n = is an absorbing state in the sense that the 
system can never leave it. According to Eq. I|ll|) any 
Hamiltonian function describing a system with empty ab- 
sorbing state must satisfy 

H R ( P ,q)\ q=o = (14) 

in addition to Eq. I|13|) . 

Before considering the full "quantum" problem, 
Eq. (|ll)p. let us analyze the corresponding classical dy- 
namics. The classical equations of the motion are 



d 

dtq = -^H R (p 7 q) , 
d 

dtp = - —H R (p,q). 



(15) 
(16) 



Due to Eqs. fTJ| and 113fl . p = 1 is always one of the 
solutions of Eq. I|16l) . Substituting p = 1 into Eq. (TSJ, 
one obtains [for the Hamiltonian (|ll|l ] 



d t q = jjy ( m ~ k ) q k ■ 



(17) 
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FIG. 1: (Color online) Phase portrait of the binary annihi- 



lation system, 2A 



The corresponding classical Hamil- 



tonian is given by Hn(p,q) = w(l — p )q ■ Solid black lines 
show zero-energy trajectories: generic lines p — 1 and double 
degenerate q — and the "accidental" line p = — 1. Dashed 
colored curves indicate trajectories with nonzero energy. The 
arrows show the evolution direction. 



This is nothing but the mean-field rate equation for the 
average particle number (n(t)}, neglecting all fluctuation 
effects. Therefore one may identify the variable q as the 
reaction "coordinate" (in fact this notation is not precise, 
since it is true only along the line p = 1). 

In order to proceed with the classical problem, 
Eqs. i|15|) and l|16fl . beyond the reaction rate approxima- 
tion (j> = 1), it is convenient to consider the phase space 
of the system. The classical equations of motion <|15[) 
and lfR)|) conserve energy. As a result, the phase-space 
(p, q) evolution of the system takes place along the tra- 
jectories given by Hji(p,q) — const, where the constant 
energy is determined from the initial conditions. Among 
all possible trajectories the ones with Hr = play a 
special role. For one thing, the evolution prescribed by 
the rate equation (p = 1) takes place along one of such 
lines [cf. Eq. More importantly, the trajectories 

with Hji = may intersect each other. Indeed, the two 
zero-energy lines, guaranteed by Eqs. and l|14|) . i.e., 
p = 1 and q = O-intersect in the point (1,0). Therefore 
the set of intersecting zero-energy curves plays the role 
of separatrix; i.e., it divides the entire phase space on 
the isolated sectors. All other trajectories cannot inter- 
sect the zero-energy ones and are confined to one of the 
sectors. The web of the zero-energy curves uniquely de- 
termines the topology of the phase portrait. An example 
of such a construction is given in Fig. ^ 

Going back to the full "quantum" problem, one may 
write a formal solution of Eq. (|10|l as 



G (Pf^f)^ J U(pf,tf,po,t )G (p )dp , (18) 

where the Green function U(pf,tf,po,to) is given by the 
T exponent: T exp{H^t}. Dividing the time interval 
[io, tf] into N — > oo steps and introducing the resolution 
of unity at each one of them, one obtains the Feynman 
representation 



J7(p/,*/;po,*o) = Vp(t)Vq(t) 



S\p,q] 



(19) 



with the Hamiltonian action 



S\p,q] = / dt\pd t q-H R (p,q)] . (20) 

J to 

To combine reaction kinetics with the random walk 
on a lattice, one needs to modify the master equation. 
The corresponding generating function becomes a func- 
tion of many variables pi, where the index i enumerates 
the lattice sites. One may also introduce the conjugated 
variables q\ = d/dpi. The resulting Hamiltonian takes 
the form 

where D is a hopping probability per unit time and the 
sum in the first term on the RHS runs over nearest neigh- 
bors Taking the continuum limitSiSa and introduc- 
ing the pair of canonically conjugated fields p(x, t) and 
q(x,t), one arrives at the quantum field theory with the 
Hamiltonian action © . The diffusion constant in Eq. © 
is given by D = Da 2 , where a is the lattice constant. 

Unless the system is very close to extinction, the func- 
tional integral in Eq. (|19|l may be evaluated in the saddle 
point approximation. In such a case the Green func- 
tion is given by the exponentiated action of a classical 
trajectory, satisfying the proper boundary conditions 69 
[much the same way as the minima of the potential V(q), 
Eq. J2J), dominate the partition function away from an 
equilibrium phase transition] . A possible phase transition 
may be associated with a qualitative change in the behav- 
ior of the phase-space trajectories (cf. the rearrangement 
of minima of the potential upon an equilibrium transi- 
tion). In other words, phase transitions lead to a change 
of topology of the phase-space portrait. Since the latter 
is determined by the set of the zero-energy lines, it is 
the rearrangement of this set, upon variation of a con- 
trol parameter, which must be associated with the phase 
transition. 

If a system is close enough to a phase transition (or 
extinction) , the saddle point approximation may lose its 
validity (below critical dimensionality d c ). One can then 
employ the RG technique to progressively integrate out 
the small scale fluctuations. Upon such a procedure the 
constants and even the functional form of the reaction 
Hamiltonian flow. However, one needs to follow the 
topology of the phase space, rather than a specific form 
of the trajectories. Around the transition the topology 
may be fully encoded in a relatively simple polynomial, 
which in turn provides a full characterization of the tran- 
sition (at least for small e = d c — d). Considering all 
distinct topologies, stable upon RG transformations, one 
may classify the possible phase transitions. 

We turn now to an illustration of these ideas with spe- 
cific examples. 
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FIG. 2: (Color online) Phase portraits of the system exhibit- 
ing the first-order phase transition. Thick solid lines repre- 
sent the zero-energy trajectories, (a)-(c) the set, Eq. I22H : (a) 
extinction phase, (b) transition point, and (c) active phase, 
(d) Reaction of the type of Eq. 12 It with k = 2, n = 3, and 
i = 4 in the active phase. 



III. MODELS WITH FIRST-ORDER 
TRANSITIONS 

Consider a set of reactions, given by 

A^>0, 2A A 3A, 3A^2A. (22) 

The corresponding reaction Hamiltonian, according to 
Eq. may be written as 

H R = (a~|A+~pV) (l-p)q. (23) 

There are three zero-energy curves: the two generic lines 
p = 1 and 9 = 0, following from Eqs. (|13[) and l|14l) . and 
the additional "accidental" curve p 2 = 2\/(fiq — aq 2 /3). 
The phase portrait of the system is depicted in Figs.^a)- 
Htc). Its topology is qualitatively different for fi < fi c = 
y/8\<r/3 and fi > fi c . In the former case, Fig. EJa), the 
"accidental" curve does not intersect the mean-field line 
p = 1 and the flow along the latter is directed towards 
q = (extinction). For fi = ji c the two curves touch 
each other, Fig. Hfb), creating a stationary point with a 
finite concentration q = 3/j/(2<t). At a larger creation 
rate /x > there is a stable point with concentration 
(n), Fig. He). 

This is the first-order transition scenario. Indeed, the 
stable concentration experiences a discontinuous jump 
from zero at fj, < ji c to some finite value at fj, > fj, c . 
It may be generalized to any reaction set of the form 

kA— >(Je - I) A , nA^(n + m)A, iA—>jA, (24) 



where the first two reactions satisfy the condition k < n 
and the last reaction is necessary to ensure a finite parti- 
cle density by having i > n,j (for hard-core "fermionic" 
models this restriction is intrinsic and the last reaction 
may be omitted). The corresponding reaction Hamilto- 
nian may be written as 

Hr = [hxijp) - h 2 (p)q n - k + h 3 {p)q l - k ] (1 -p) q k , (25) 

where the functions hi,h,2, and ft. 3 are positive in the 
interval p € [0, 1]. An example of the phase portrait in 
the active phase and k — 2 is depicted in Fig.J^d). 

Strictly speaking, the condition /i = // c [Fig. H£b)] 
signifies the appearance of the metastable state (unless 
D — » oo). Such a metastable state becomes stable at 
some larger creation rate jl c {D) > /i c . Therefore, in 
terms of the bare system's parameters the actual first- 
order transition takes place when the phase portrait has 
a form of Fig. Etc). Alternatively, one may imagine in- 
tegrating out spatial fluctuations (governed by the diffu- 
sion constant D) and plotting the phase portrait of the 
effective zero-dimensional (d = 0) system in terms of the 
renormalized parameters. It is in this latter sense that 
the transition point is depicted by Fig.EJb). 

An interesting limiting cas o 51 i 52 of the first-order tran- 
sition is the reaction set 

kA^(k-l)A, kA-^{k + m)A, (26) 

described by the Hamiltonian 

H R = h(p)(l-p)q k , (27) 

where h(p) is a polynomial the degree k + m — 1. The 
zero-energy lines are given by p — 1, the k times de- 
generate line q = along with the lines p = pi, where 
Pi are roots of the polynomial h(pi) — 0, Fig. |3J It is 
easy to check that for \x c — one of the roots of the 
polynomial is pi = 1. In this situation the p — 1 zero- 
energy line is doubly degenerate, corresponding to the 
first-order transition. Upon such a transition the p = 1 
and p = p\ lines interchange their relative positions and 
the extinction behavior turns into the unlimited prolifer- 
ation one. 

Notice that adding a proliferation-restricting reaction 
iA — > jA, with i > k, to the set (I26|l qualitatively 
changes its phase portrait. In general, such a modifi- 
cation leads to a continuous phase transition [see, e.g., 
Eqs. I|28|l and (|38|l below]. An analogous transformation 
of the first-order to the continuous transition takes place 
for restricted ("fermionic") reaction rules22i31. This is 
consistent with the conjecture that the "fermionic" con- 
dition implies an effective annihilation process of the or- 
der of k + 1; see, e.g., Refs. and references therein. 

The phase portrait description given above has the sta- 
tus of the mean-field consideration only (much the same 
way as the finding of minima of the effective potential, 
Eq. (J2J, in the equilibrium theory). It does not exclude 
the possibility that below some lower critical dimension 
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FIG. 3: (Color online) Phase portraits of the system, Eq. 1261 . 
exhibiting the first-order phase transition. Thick solid lines 
represent the zero-energy trajectories, (a)-(c) Reaction set 
with k = I = m = 1: (a) extinction phase, (b) transition 
point, and (c) unlimited proliferation phase, (d) Reaction set 
with k = I — 2 and m = 1 in the extinction phase. 



FIG. 4: (Color online) Phase portrait of the DP system in 
the active phase. Thick solid lines represent zero-energy tra- 
jectories which divide the phase space into a number of dis- 
connected regions. Point B = {1, 2(u — A)/cr} represents the 
active mean-field point. The system is brought to the phase 
transitions if points A, B, and C coalesce. 



d' c the fluctuations may drastically change the system's 
behavior. In particular, a transition, being first-order in 
the mean-field scenario, may be turned to a continuous 
one by fluctuation effects. This may explain the appar- 
ent continuous transition seen in d — 1 systems, similar 
to those given by Eq. (|24f! 49 - 71 ' 72 . It was also observed 
in some fermionic reaction schemes that the first-order 
transition observed at large D turns out to be continu- 
ous for D < D,m2i21. The subject definitely needs more 
investigation both analytical and numerical. 

IV. DIRECTED PERCOLATION AND ITS 
GENERALIZATIONS 

A. DP models 

Consider a reaction set which includes death, branch- 
ing, and coalescence reactions: 

A^>0, A^2i, 2A^A. (28) 

The corresponding reaction Hamiltonian takes the form 

Hr = X(l-p)q + p{p 2 -p)q+^(p~p 2 )q 2 

= (a-mp+I^?) {i-p)q- (29) 

The phase portrait of the DP system is depicted in 
Fig. ^ The lines of zero energy are generic p = 1 and 
q = trajectories along with the "accidental" trajectory 
q = 2((ip — X)/ap. According to the mean-field analysis 
[classical equations (|15|) and (|16fl with p = 1] , there is an 
active phase with the average density 

(n) = 2 ^— - (30) 
a 

for p > A. The active phase corresponds to point B 
in Fig. 2| The system may be brought to extinction by 



driving the control parameter m = /i — A through zero. 
Therefore the p = A point corresponds to the continuous 
phase transition. The transition is represented by the 
phase portrait with the three zero-energy trajectories in- 
tersecting at the single point (1,0). According to Eq. l|3Tj|) 
the mean- field order-parameter exponent is = 1. The 
other mean-field exponents^ are v± — 1/2, vu = 1. 

To go beyond the mean-field picture one needs to in- 
vestigate the immediate vicinity of the phase transition. 
To focus on this regime it is convenient to shift the mo- 
mentum variable 

p-l->p. (31) 

Moreover, close to the transition the phase portrait and 
thus the Hamiltonian may be modeled by the three in- 
tersecting straight lines, Fig. Efa)- EJc). This way one 
arrives at the model Hamiltonian, applicable close to the 
DP transition: 

Hr = (to + up — vq) p q . (32) 

The bare values of the constants are given by m = (J, — A, 
u = fi, and v = a/2. The corresponding action, Eq. 
acquires the formiSii 

S = Jdt d d x [p(dt — -DV 2 ) q — rap q — up 2 q + vp q 2 ~\ , 

which may be recognized as a Reggeon field theory 
actionZ^. 

There are many other reaction sets, in addition to 
Eq. H28|) . with the same "triangular" topology of the 
phase portrait. Some of the examples are A — > ; A — > 
(to + I) A ; 2A -* A, m > 1 , or A ->• 2A ; 2A -> 
(see also Sec. lIVCTl : etc. In the vicinity of the phase tran- 
sition they all exhibit the topology of the phase portrait 
depicted in Figs. EJa)- E{c). Therefore they all may be 
described by the model Hamiltonian (|32|l . Accordingly 
they all belong to the same DP universality class. 



FIG. 5: Generic phase portrait of DP models in the vicinity 
of the phase transition [after the shift, Eq. (|31|t]. (a) Active 
phase, m > 0; (b) transition point, m — 0; (c) extinction 
phase, m < 0. The plus and minus signs show the sign of 
the Hamiltonian in each sector, (d) The one-loop diagram 
renormalizing u vertex (vertices m and v are renormalized in 
a similar way). 

Naive scaling dimensions of the action (|33[1 are z = 2, 
while [p] + [q] = d. Since one expects^ both vertices u 
and v to have the same scaling dimensions, one finds [p] = 
[q] = d/2. As a result, the bare scaling dimensions of the 
vertices are [m] = 2, while [u] = [v] = 2 — d/2. Therefore 
below the critical dimension d c = 4 the nonlinear vertices 
u and v are relevant and the mean-field treatment is 
expected to fail. 

The one-loop corrections to the naive scaling are given 
by the triangular diagrams, like the one depicted in 
Fig. Efd). Such diagrams are logarithmically divergent 
in d = 4, as expected. Straightforward calculations^ 
(see Ref. ^] for the two-loop approximation) lead to 
the following set of RG equations: 

dim = (2-Suv)m, (34) 
diu = (e/2-2Suv)u, (35) 
div = (e/2-2Suv)v, (36) 

where e = 4 — d and the differentiation is over the log- 
arithm of the scaling factor. We have introduced the 
factor 

s-a 1 f A ddk A '~ 4 r,7i 

which may be absorbed in the proper redefinition of the 
running constants. 

According to Eqs. (|3~5j) and (|3*dT) . di(u/v) = 0, mean- 
ing that the slop of the "accidental" zero-energy line, 
q = (m + up)/v, remains intact upon the renormal- 
ization procedure. In fact, this statement is exact be- 
cause of the symmetrjiS. As a result, the overall topol- 
ogy of the phase portrait is preserved by the RG. For 
d < 4 the RG equations (|34|l - 136 (1 predict a nontriv- 
ial fixed point given by m* = 0, u* = \/ efi/(ASa), 
and v* = y/ e<r/(4Sn). Substituting these values into 
Eq. I|32|) . one finds the fixed-point reaction Hamiltonian 



FIG. 6: (Color online) Phase portrait of the 2CPD system in 
the active phase (cf. DP, Fig. |1J. The zero-energy line q = 
is doubly degenerate and is depicted by the double line. At 
the transition points A, B, and C coalesce. 



Hft(p, q) corresponding to the DP universality class phase 
transitions. Its phase portrait is depicted in Fig. E£b). 
Linearizing the RG equations l|34|) - (|36|l near the fixed 
point, one finds dim = v^m with the critical exponent 
v± = (2 - e/4)" 1 w 1/2 + e/16. The other critical ex- 
ponents may be deduced in the standard wayii, e.g., 
/3« l-e/6. 



B. fc-particle contact processes 

As mentioned in the Introduction, the DP universality 
class is extremely robust. This fact is due to the stability 
of the "triangular" topology of the phase portrait near 
the transition. One may try to change this topology by, 
say, requiring four or more zero-energy trajectories to 
intersect in the same point. It is clear, however, that 
in general one must fine-tune more than one parameter 
to reach such a scenario. Even if bare reaction rates 
are specially chosen to let it happen, the tuning is not 
expected to survive upon RG integration of fluctuations. 
Therefore crossing of more than three lines is possible 
only at a multicritical transition point. The only way to 
go beyond the DP is if a different topology of the phase 
portrait is imposed by an additional symmetry. 

In this section we discuss a class of models where the 
minimum number of particles needed to initiate any reac- 
tion is restricted to be k > 1. According to Eq. (fTT|) all 
terms in the corresponding reaction Hamiltonian must 
contain q k or higher power. In terms of the phase por- 
trait it means that the generic q = zero-energy trajec- 
tory is k times degenerate. To emphasize the difference 
with nondegcncratc lines, we denote such fc-degenerate 
trajectories by k closely spaced parallel lines. An impor- 
tant fact is that the degeneracy is preserved by the RG 
transformations. Indeed, the fluctuations cannot initiate 
a reaction with fewer than k incoming particles, if it is 
not in the original reaction set. We denote such models 
as fc-particle contact processes with diffusion (kCPD). 
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To be specific, let us consider the case of k = 2, which 
is represented by, e.g., the following set of reactions: 



2A 



2 A A 2, A ■ 3A-+2A. 



(38) 



The corresponding reaction Hamiltonian takes the form: 



A 
2' 



II n = ~ ( i - p 2 ) g 2 + f 3 - ) + ^ (p 2 - p 3 ) g 3 



u 
6< 



\ (\(l+p)-VP 2 + ^p 2 q) (l- P )q 2 



(39) 



The phase portrait of the 2CPD system is depicted 
in Fig. The lines of zero energy are generic p = 1 
and double-degenerate q = trajectories. The "acciden- 
tal" trajectory is given by q = 3 \ji,p 2 — A(l +pj\ /(o~p 2 ) ■ 
There is an active phase with average density (n) — 
(3/U — 6A)/er for fi > 2A. It corresponds to point B in 
Fig. The system may be driven to extinction by tun- 
ing the control parameter to = /i/2 — A to zero. Therefore 
the n/2 — A point corresponds to the continuous phase 
transition. At the transition point the four (we count the 
q = line twice) zero-energy lines are intersecting in the 
single point (1,0). 

Focusing on the transition region, one may shift the 
momentum variable p — 1 — > p and model the zero-energy 
trajectories by straight lines, Fig. 0a). The resulting 
model Hamiltonian, applicable close to the transition, is 



Hff = (to + up — vq) p q 2 



(40) 



The bare values of the constants are given by m = n/2 — 
A, u = fj, — A/2, and v = a/6. Since near the transition 
to ss and thus /i ss 2 A > 0, one finds for the bare value 
m w 3A/2 > 0. Apart from Eq. (|38|l . there are infinitely 
many other reaction sets, which have the same topology 
of the phase portraits. Therefore the phase transition 
of these other models is described by the same model 
Hamiltonian 1 (301 . Examples include 2 A — > A , 2 A — » 
AA, AA -> 0, etc. 

In an analogous way, one may show that the phase 
portrait of a generic kCPD process, such as, e.g., kA — > 
, fe^4 — > (fe + 1)A , (fe + 1)A — > fc/1, contains a triangle of 
fc-times-degenerate q = line, p = line [after the shift, 
Eq. H31|) ]. and the "accidental" q = (m + up)/v diagonal 
line; see Fig.[7fb). Thus it may be described by a model 
Hamiltonian of the form 



Hr = (to + up — vq) p q k 



(41) 



where to is the control parameter of the transition and u 
and v are positive constants. 

To find scaling exponents near the transition one as- 
signs bare dimensions z — 2 and [p] + [q] = d. There are 
no perturbative corrections to the propagator (for k > 2), 
and thus one does not expect these exponents to be 
changed in the e expansion. Since both u and v vertices 
have to be equally (ir)relevant on the mean-field level, 
one has to choose the bare dimensions as [p] = [q] = d/2. 
As a result, the bare scaling dimensions of the vertices 




7a 



(b) 



FIG. 7: Generic phase portraits of (a) 2CPD models and (b) 
3CPD models in the active phase. 



are [m] = 2 — (k — l)d/2, and [it] = [v]=2- fed/2. Ac- 
cording to this scaling analysis the critical dimension is 
expected to be d c = A/k. The corresponding mean-field 
exponents of the kCPD transitions at d > d c are (3=1 
and u±_ = k/2. Only the k = 2 and k — 3 processes 
are expected to exhibit nontrivial behavior in the phys- 
ically relevant dimensionsSiSS. We shall analyze these 
two cases below. 

The 2CPD transition is described by the action © 
with the reaction Hamiltonian l|40|) (see, e.g., Ref. y|). 
Its critical dimension is d c = 2. The one-loop renormal- 
ization is given by two- vertex loops, which are logarith- 
mically divergent in d — 2. One arrives^ at the following 
set of RG equations: 



dim = (1 + e/2 + Su)m, 
diu = (e + Su) u , 
div = (e + 3Su)v, 



where e = 2 — d and 



s = a 



At: 



d d k 

-i 15k 2 



A 



d-2 



AttD 



(42) 
(43) 
(44) 



(45) 



Notice that the sign of the perturbative corrections in 
Eqs. (|42|l - l|44|l is opposite to that in the corresponding DP 
RG equations l(3~4l - ((3l)|l . As a result, the weak-coupling 
fixed point appears to be absolutely unstable for d < d c 
(for positive initial u). Its region of stability for d > d c 
is finite and shrinking as d — > d c from above^i. Solving 
Eq. one finds that the coupling constant u diverges 
once the RG reaches a certain spatial scale £ -1 = Ae~', 
where, in d < 2, 



l/(2-d) 



(46) 



(here uq is an initial value of u). In d — 2 one finds 
that £ = A -1 exp(47r_D/uo)- This indicates that some 
new physics shows up at the scale £. In Appendix lAl we 
suggest that the system may develop anomalous aver- 
ages, similar to those in BCS theory. The corresponding 
"coherence length" appears to be exactly £. (See also 
Refs. [8u67j for further discussion.) 

A similar situation is encountered for k = 3. The crit- 
ical dimension is d c = 4/3. The RG analysis of Eq. I|41|) 
with k = 3 shows that the only vertex which acquires per- 
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turbative corrections is v. The corresponding RG equa- 
tions are 



dim = (2 — d) m , 
dm = (2-3d/2)«, 
div = (2-3d/2)v-3Smu. 



(47) 
(48) 
(49) 



Once again, the coupling constant u grows indefinitely 
and the e expansion fails to predict critical exponents. 








active 



empty 



C. Critical point and mean— field transitions 

In the discussions above we have avoided considering 
an exact location of the critical point in the parameter 
space of the problem. Instead, we have concentrated on 
the universal properties of the transition itself. The only 
thing we had to assume is that the three intersection 
points A, B, and C in Figs. |S] and can be made to 
collapse. In other words, the parameter m in Eqs. 1-iL'l) 
or l|41|) can be made to be zero. In terms of the bare 
parameters of the model [e.g., Eq. I|28(l ] this corresponds 
to /i = fcA. However, the actual transition may take 
place away from the (i = k\ point. Indeed, it is the 
renormalized parameter m that vanishes at criticality, not 
the bare one. 

The particular case of bare A = 0, e.g., 



kA 



(k + l)A, (k + l)A 



(50) 



has attracted a lot of attentiona*iaiAi. The naive mean- 
field expectation is that such a system is always in its 
active phase, unless [i — 0. Yet it is clear that any finite 
size system is bound to end up in its empty absorbing 
state after long enough time. [The latter may be esti- 
mated as an exponentiated classical action accumulated 
along the trajectory leading from the active mean-field 
state (point B in Figs. 0] and HJJl to the empty absorb- 
ing state (point Therefore, if the diffusion time 
between the sites is longer than the respective extinc- 
tion time, the system evolves towards extinction even at 
fi > 0. As a result, there is a nontrivial phase separation 
boundary in the parameter space (D,n) of the A = 
systemic, Fig. |SJ One may interpret the appearance of 
this nontrivial line as the generation of the effective anni- 
hilation rate A e //(-D, /j.) > 0. For example, two successive 
processes A — > 2A and 2A — > result in an A — > ef- 
fective rate. Consequently the phase transition along the 
nontrivial critical line is represented by the evolution of 
the phase portraits depicted in Figs. and and belongs 
to the kCPD universality classes. 

The crucial observation made by Cardy and Tauber— is 
that at and below a certain lower critical dimension d' c the 
separation boundary extends to an arbitrarily large dif- 
fusion constant D, Fig.|SJa). This fact may be associated 
with a divergent perturbative correction to the coupling 
constant m. In the two most interesting cases DP and 
2CPD such a correction is given by the two-vertex loop, 



FIG. 8: Schematic phase boundary in the parameter space 
(D, /x) of the A = reaction set (I5UH . (a) d < d' c , the transition 
along the entire line is of kCPD class, (b) d > d' c , to the left 
of the tricritical point T the transition is of kCPD class, while 
to the right it is /i = mean-field transition. 




FIG. 9: (Color online) Phase portrait evolution upon /i = 
mean-field transitions: (a) k = 1 and fj, > 0, decreasing 
the branching rate fi is schematically shown by the dashed 
arrows, (b) k = 1 and \i = 0, the "accidental" line becomes 
horizontal and the local topology of the phase portrait is that 
of the binary annihilation reaction; cf. Fig. (c) k = 2 and 
(i > 0, (d) k = 2 and /i = 0, the topology is that of the 3-plet 
annihilation reaction, 3^4 — » 0. 



Eq. 1)45(1. As a result, the corresponding lower critical 
dimension is d' c — 2. Above the lower critical dimen- 
sion, d > d' c , there is a tricritical point T = (D c , 0), such 
that for D > D c the transition occur only at /i = 0^, 
Fig-IHtb). 

The evolution of the phase portrait corresponding to 
such a fj, = transition is depicted in Fig. |5| (cf. Figs. 
El or 01 . At /i = the only remaining reaction is n > 2 
particle annihilation, (k + 1)A — > 0, with the critical 
dimensionality d c = 2j(k — 1) < 2. Thus, the critical 
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dimensionality of the /i = transitions is smaller than 
d' c — the minimum dimensionality above which the fi = 
transitions may take place. As a result, all such transi- 
tions occur above their respective critical dimensionality 
and therefore exhibit mean-field behavior. The generic 
scheme for a mean-field fi = transition may be written 
asffi 

kA^{k + l)A, nA-^(n-m)A, (51) 
where k < n. The corresponding rate equation is 

d t q = r t -Lq k -a-q n . (52) 
k\ n\ 

The stationary density is proportional (n) oc 
(/Lt/cr) 1 / ( -™ _fc - ) , leading to the mean-field critical ex- 
ponent j3 = l/(n — k). The other mean-field exponents 
can be also deduced in a straightforward waj^S. 

There are two important exceptions where the fj. = 
transitions in the reaction set Eq. (|51|l may exhibit non- 
mean-field behavior. These are the cases where genera- 
tion of the effective A e ff(D,(i) is prohibited by certain 
symmetries. This may lead to a scenario with d' c < d Cl 
thus leaving a window d' c < d < d c for non-mean-field 
fi = transitions. There are two such symmetries: (i) 
parity conservation (both I and m are even), discussed 
in Sec. [V] and (ii) reversal symmetry (k + I = n and 
n — m — k), discussed in Sec. IVII 

V. PARITY CONSERVING MODELS 

A. PC model 

As was mentioned above the parity conservation may 
lead to a new nontrivial universality class. Consider the 
simplest possible PC reaction set 

A A 3/1 , 2,4 A 0. (53) 

The corresponding reaction Hamiltonian takes the form 

H R = {up-vq){p 2 -l)q, (54) 

where the bare values of the constants are given by u — 
fi and v = a/2. The corresponding action is invariant 
under the following transformations 5 : 

P^-P, ?-►-<?> (55) 

which may be traced back to the conservation of par- 
ity. As a result, the phase portrait, Fig. ^3 possesses 
the reflection point at the origin. This symmetry is pre- 
served upon RG transformations. Therefore the "acci- 
dental" zero-energy line [q = up/v, according to Eq. I)54p] 
is bound to be an odd function and thus pass through the 
origin. Its shape, however, may change in the process of 
renormalization. Consequently the phase transition can- 
not be described by the coalescence of three points A, B, 
and C and its nature is different from the DP class. 
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FIG. 10: (Color online) Phase portrait of the parity conserv- 
ing model, Eq. (I53H . Notice the reflection symmetry around 
the origin. 



According to the mean- field equation dtq = 2fiq — aq 2 , 
the model is always in the active phase with the number 
of particles (n) = 2/i/a = u/v. The only way to drive 
the mean-field dynamics towards extinction is to send 
u = (J, — ► 0. In other words, the critical point is u c = 0. 
One may discuss, though, the scaling of particle density 
with u — u c — u and this way define the "magnetization" 
exponent (3: 

(n) - u , (56) 

where the mean-field value of the exponent is /3 = 1. 

Turning to the fluctuations, one notices that it is not 
possible to perform the shift of momentum, Eq. I|31|l . and 
focus on the immediate vicinity of the (1, 0) point on the 
phase plane. Because of the symmetry (|55|) one has to 
keep the entire interval p G [—1,1] under consideration; 
see Fig.^3 Therefore one must choose the scaling dimen- 
sion [p] = 0. Since the bare scaling requires [p] + [q] = d, 
the naive scaling dimension of q is [q] = d. As a result, 
one finds (since z = 2) [u] = 2 and [v] = 2 — d and the 
critical dimension is d c = 2. One may worry that since 
[p] = 0, it is not possible to restrict the consideration 
to the low-order polynomial in p. Instead, one has to 
keep all the powers of p, resorting to the functional RG. 
We perform such a procedure in Appendix [B] and show 
that it actually justifies the use of the truncated reaction 
Hamiltonian l|54l) . 

The one-loop RG calculation, utilizing the two-vertex 
loop (which is logarithmically divergent in d = 2), yields 
the following RG equations^: 

diu = (2-3Sv)u, (57) 
d t v = {e-Sv)v, (58) 

where e = 2 — d and S is given by Eq. I|45|) . For e > 
there is a nontrivial stable fixed point v* = e/S. In the 
vicinity of this fixed point the relevant parameter u scales 
as diu = (2 — 3e) u and thus its new scaling dimension is 

[u] = 2 - 3e = id - 4 . (59) 



This leads to the non-mean-field exponent (3 given by 

= d/[u] « 1 + e. (60) 

The fact that (3 > (3 MF = 1 means that the actual den- 
sity in d < 2 is less than the mean-field prediction. The 
fluctuations drive the system closer to extinction. Cardy 
and Tiiuber— suggested that for d < d ' c « 4/3 there is an 
extinct phase at finite \x and the transition to the active 
phase takes place at some /i c > and estimated the cor- 
responding critical exponents. Recently, the PC model 
was reexamined employing the so-called nonperturbative 
renormalization group method by Canet et tdmi±-*. They 
confirmed the existence of the /j, > transition for the 
PC model below d' c = 4/3 and estimated the exponents 
with fair accuracy for d = 1 . 

B. Generalized PC models 

One may invent other models conserving parity. For 
example, Cardy and Tauber— considered the class of par- 
ity conserving models 2A — > 0, and A — > (m + l)A 
with even m. The corresponding reaction Hamiltonian 
is H R = [uph m (p) - vq]{p 2 - l)q, where h m (p) = 
(p m — l)/(p 2 — 1) is an even polynomial. Its phase portrait 
is topologically identical to Fig. El Thus one expects this 
reaction set to belong to the same universality class as 
the PC model. This was indeed the conclusion of Ref. 

To find a different topology of the phase portrait and 
therefore a new universality class one needs to impose 
an additional symmetry. Following Sec. IIVBI we shall 
consider parity conserving reactions with minimal num- 
ber fc > 1 needed to initiate all reactions. For example, 
consider a parity conserving set of reactions with even fc: 

kA^(k + 2)A, (k + l)A A A, (61) 

The reaction Hamiltonian is 

H R = [up"- 1 - v h k (p) q] p(p 2 - 1) q k , (62) 

where hk(p) is an even polynomial of the degree fc — 2, 
u = /i/fc!, and v — a/(k + 1)!. The corresponding phase 
portrait is depicted in Fig. Illf a). For an odd k a repre- 
sentative set of reactions is 

kA^(k + 2)A, (fc + l)A^>0. (63) 

with the reaction Hamiltonian 

H R = [up k - v h k+1 (p) q) (p 2 - 1) q k , (64) 

where hk+i(p) is an even polynomial of degree k — 1. The 
corresponding phase portrait is depicted in Fig. ITTT bh 
These phase portraits are topologically stable upon RG 
transformations and thus represent a set of distinct uni- 
versality classes. We call them fc-particle parity conserv- 
ing (kPC) classes. 
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FIG. 11: Phase diagrams of kPC models, (a) 2PC model, the 
inset shows the diagram leading to the logarithmic corrections 
in d = 1. (b) 3PC model, (c) 2PC model with all reactions 
including only an even number of particles. This topology is 
not stable and evolves towards (a), (d) 2PC with 2 A — > 0, the 
corresponding topology is essentially equivalent to the 2CPD 
model, Fig. El 

Assigning the scaling dimensions as in the PC model, 
[p] = 0, [q] = d, and z = 2, one finds [u] — 2 — (fc — l)d 
and [v] = 2 — kd. At the critical dimension v turns out to 
be relevant and thus d c = 2/fc. Therefore at any physical 
dimension the kPC behavior is described by the mean- 
field. The only exception is the 2PC model, Fig. IllT ah 
which acquires logarithmic corrections to the mean-field 
scaling at d = 1. The renormalization is due to the two- 
loop diagrams built with the help of the p 3 q 3 vertex; see 
the inset in Fig. Illf ah 

Other attempts to generalize the PC universality class 
appear to be unstable against RG transformations. For 
example, consider a parity conserving set which contains 
only an even number of reagents: 2 A — > A A and AA — > 
0. The corresponding reaction Hamiltonian is given by 
Hp = [up 2 — v(l +p 2 )q 2 ](p 2 — l)g 2 . Its phase portrait is 
depicted in Fig. Illf ch In addition to the PC symmetry, 
Eq. H55|) , the Hamiltonian and the phase portrait possess 

q-> -q (65) 

symmetry. However, this is not the symmetry of the full 
action, Eq. J3J. Therefore this symmetry is not stable 
against the RG transformations. Indeed, e.g., using three 
vertexes p 2 q 2 , one may generate a p 3 q 3 vertex which vio- 
lates the symmetry (|65|) . As a result the system belongs 
to the 2PC class and its phase portrait drifts towards 
Fig. Eta). 

One may add a competing annihilation reaction to the 
kA — > (fc + 2) A process of kPC for fc > 2, such as kA — > 
(fc — 2) A. For example, 

2AA4A, 3A^A, 2A^0. (66) 

Because of the competition, one expects the absorbing 
state transition to happen at m — 2/i— A = 0. The Hamil- 
tonian is Hr = (m — u + up 2 — vpq){p 2 — l)q 2 ; it obeys 
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the PC symmetry l|55fl. The phase portrait is plotted in 
Fig. Illf dh One notices that in the vicinity of the tran- 
sition the local topology is indistinguishable from Figs. El 
and Ufa). Therefore the transition belongs to the same 
universality class as 2CPD (see Sec. IIV Bl) . This fact 
was already noticed in numerical simulations^*^, but, to 
the best of our knowledge, remained unexplained. The 
identity of the phase portrait's topologies in the vicinity 
of transition immediately explains the universality. One 
can show that the other kPC processes with competition 
in the fc-particle channel belong to the same universality 
classes as the corresponding kCPD models. The example 
of the model, Eq. ljf)6")l . shows that parity conservation is 
not a crucial feature to discriminate the PC class (see also 
Refs. It is rather the topology of the phase por- 

trait in the vicinity of the transition that discriminates 
various universality classes. 



VI. REVERSIBLE REACTIONS 
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FIG. 12: (a)-(c) Phase portraits of reversible reactions, (a) 
A <-> 2A and (b) 2A <-> 3A. (c) Parity conserving A <-> 3A. 
(d) A topology with q — > —q symmetry unstable against RG 
transformations; see Sec. I VIll 



All the phase transitions, discussed above, are associ- 
ated with a deformation and rearrangement of a charac- 
teristic triangular structure (possibly with one degener- 
ate side) on the phase plane. There is one more possibil- 
ity for a stable transition topology which is a rectangular 
structure. All models exhibiting rectangular topology 
consist of a single reaction which is allowed to go in both 
directions with different rates. 

Consider, for example, the reversible reaction [same as 
DP, Eq. QgJ, but without A -> 0] 

AA2i, 2A^A. (67) 
The corresponding reaction Hamiltonian is 

= (p 2 -p)(uq-vq 2 ), (68) 

where u = /i and v = a/2. The phase portrait is depicted 
in Fig. I12f a) and has characteristic rectangular shape 
comprised by the generic lines p — 1 and q = along with 
the two "accidental" ones p = and q = u/v. The mean- 
field predicts the average density to be (n) = u/v. The 
control parameter is u — \i with the critical value u c = 0. 
Consequently the mean-field "magnetization" exponent 
is j3 = 1, in agreement with the general schemed; cf. 
Eqs. (|51|l and 152(1. One may ask if the mean- field be- 
havior can be modified by the fluctuations. 

To answer this question, one notices that the phase 
portrait topology is stable against renormalization. I.e., 
no terms violating the rectangular structure are gener- 
ated. It may be checked either by considering possible 
diagrams or realizing that the action possesses the sym- 
metry 

v u 
p->-q, q->-p, t^-t. (69) 
u v 



To keep the entire interval p £ [0, 1] unchanged upon 
rescaling^ 4 , one has to choose the scaling dimensions as 
[p] = and [q] = d, then [u] = 2 and [v] = 2 — d. From 
here one concludes that the critical dimension is d c = 2. 
The RG equations are 

d t u = {2-Sv)u, (70) 
div = {e-Sv)v, (71) 

where e = 2 — d and S is given by Eq. I|45|). For e > 
there is a stable fixed point at v* — e/S. In its vicinity u 
scales as diu = (2 — e) u\ thus its new scaling dimension 
is [u] = 2 — e = d (see also Appendix lB|l . As a result, the 
exponent is given by (3 = d/[u] = 1 + 0(e 2 ). At least in 
this order, it is not affected by the fluctuations. In fact, 
P = 1 is proven to be exact in d = 1, Ref. 0|- Using 
detailed balance, Ref. [84[ argues that j3 = 1 is exact in 
any dimension. 

A generic reversible reaction kA ^> mA, with m > k, is 
described by Hr — (p m ~p k )(uq k —vq m ). Its zero-energy 
lines are q — and p — 0, both k times degenerate, along 
with nondegenerate p — 1 and q = (n) = (it/v) 1 /*™ 1- *") 
For the parity conserving case (m — k even) also p = — 1 
and q = — (n) , Fig. I12f c) . By a proper rescaling of p and 
q the Hamiltonian may be brought to the symmetric sep- 
arable form Hr = —f(p)f{q)- The corresponding action 
is symmetric against p <-> q and t — » —t. Therefore the 
rectangular structure is stable in the course of renormal- 
ization. The topological structure is fully determined by 
the index k and the parity. 

Therefore, one may identify two more families of uni- 
versality classes. One is parity nonconserving, repre- 
sented by the reversible reaction 

kA^(k + l)A (72) 

(a higher number of offsprings, fc+l + 2n, does not change 
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the universality class) with the reaction Hamiltonian 

H R = (p-l)(u-vq)p k q k . (73) 

We denote it as fc-particlc reversible (kR). The action 
possesses the symmetry (|69|) . rendering stability of the 
rectangular topology. The upper critical dimension is 
d c = 2/k. An example is 2A <-> 3A [see Fig.lT2Thl]. with 
d c = l. 

The parity conserving reversible reactions 

kA <-> (k + 2) A (74) 

(a higher number of offsprings, k + 2n, does not alter the 
universality class) with the reaction Hamiltonian 

H R = (p 2 -\){u-vq 2 )p k q k . (75) 

We denote it as kRPC. The corresponding action is sym- 
metric against two symmetries, Eqs. I|55[) and Q69[) . They 
impose stability of the rectangular topology symmetric 
with respect to reflection around the origin. The critical 
dimension of the kRPC family is d c — 2/(k + 1). An 
example is A 3 A [see Fig. H2f c)]. with d c = 1. All 
other reactions of this type have d c < 1 and thus are 
fully described by the mean-field treatment (see, how- 
ever, Ref. 0). 

VII. CONCLUSIONS 

We have argued that the universality classes of phase 
transitions in reaction-diffusion models may be classified 
according to the topological structures of the correspond- 
ing phase spaces. This structure is fully encoded in the 
web of the zero-energy trajectories. The simplest and 
most stable structure is given by three mutually inter- 
secting lines. By changing a single control parameter the 
three intersection points may be made to coalesce; see 
Fig. |SJ At such a value of the control parameter the 
system undergoes a phase transition into the absorbing 
phase. The corresponding universality class is known as 
directed percolation. 

There are only a limited number of ways to organize a 
continuous phase transition, governed by a single control 
parameter, which utilizes topology different from the DP. 
We have identified four families of such transitions, which 
are stabilized by an additional symmetry or symmetries. 

(i) A generic reaction set constrained by the require- 
ment that all reactions need at least k incoming particles. 
The corresponding phase portrait is bound to have the 
q = line to be k times degenerate, Fig. [7] This prop- 
erty is robust against RG transformations. Indeed, no 
vertices with fewer than k external q "legs" can be gen- 
erated. The triangular topology with one /c-degenerate 
line, Fig. defines a family of universality classes. We 
call them kCPD's (k = 1 is DP). The scaling consid- 
erations suggest that their upper critical dimension is 
d c = 



(ii) A set of reactions which conserve parity. In this 
case the Hamiltonian and the action are invariant under 
the transformation 1)550. It dictates the reflection sym- 
metry of the corresponding phase portraits, Figs. ITU1 and 
1111 The symmetry is preserved upon renormalization. In 
addition to parity conservation one may require a min- 
imal number k of incoming particles for every reaction. 
This generates the phase portraits depicted in Fig. ITT1 
There is one universality class for every k, termed kPC, 
Figs. lllf al and II If b) (k = 1 is PC). Their upper critical 
dimension is d c = 2/k. To realize a kPC transition, the 
reaction starting from the minimal number k of particles 
must go only up, e.g., kA — > (k + 2) A. By adding down- 
going reaction, e.g., kA — > (k — 2) A (for k > 2) - the 
model is transformed into the kCPD class (despite of the 
parity conservation). 

(iii) A single reaction which is allowed to go both di- 
rections, with different rates: kA «-» [k + 1 + 2n)A. The 
corresponding reaction Hamiltonian and the action are 
symmetric under the exchange transformation, Eq. 1G9I1 . 
The phase portrait has the stable rectangular structure 
with /c-degenerate p = and q = lines, Figs. I12f a) and 
rT^T b). Upon decreasing the creation rate the rectangle 
collapses onto the interval p = [0, 1]. We call such tran- 
sitions kR. Their critical dimension is d c = 2/k (same as 
kPC, but exponents are different in d < d c ). 

(iv) A single reaction which is allowed to go both di- 
rections and conserves parity: kA <-> (k + 2n)A. The 
corresponding reaction Hamiltonian and the action are 
symmetric under the two symmetry transformations, 
Eqs. (|55|l and (|69[) . The corresponding phase portrait 
has the stable "checkered" structure, Fig. I12f . which 
collapses when sending the creation rate to zero. We de- 
note such transitions as kRPC. Their critical dimension 
is d c = 2/(k+ 1). 

Altogether we identify five nontrivial universality 
classes with d c > 1: DP, PCPD, 3CPD, PC, and 1R. 
In addition there are four marginal classes with d c = 1: 
4CPD, 2PC, 2R, and 1RPC. In sufficiently high dimen- 
sions there are other universality classes not belonging 
to the four families described above. However, they are 
bound to be of the mean-field type and do not have non- 
trivial representatives. 

We have not found any other stable, topologically dis- 
tinct structures in the phase space. Consider, for ex- 
ample, the topology depicted in Fig. IT2T d'). The phase 
portrait is symmetric under the transformation q — ► —q. 
A corresponding reaction set consists of reactions all 
starting from odd number of particles: e.g. A — > 0, 
A — > 2A, and 3A — > 0. The corresponding Hamilto- 
nian [after the shift (|31|l and neglecting irrelevant terms] 
is Hr = (to + up — vq 2 )pq. By changing the parameter 
m, e.g., by increasing the rate of annihilation A — > - 
one may bring the system to the phase transition into the 
extinction phase. Naively, such a transition is associated 
with the vertex of the parabola crossing the point (1, 0); 
see Fig. ll2f d). Such a topology is different from all those 
considered above and could represent a new universality 
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class. 

However, the q — » — q symmetric phase portrait, 
Fig. I12f d). is unstable upon renormalization. Indeed, 
combining two vertices pq 3 and p 2 q in the loop, one gen- 
erates the p q 2 term, which violates the symmetry. This 
term represents the induced reaction 2 A — > A, originat- 
ing from two reactions A — > 2A followed by 3A — ► 0. In 
other words, g — ► — g is no< a symmetry of the action and 
therefore it is not preserved by the RG. As a result, the 
initially symmetric zero-energy line p — (vq 2 — m)/u is 
shifted and deformed upon renormalization. The topol- 
ogy drifts either towards the first-order scenario, Fig. [21 
or towards the DP, Fig. [SJ To keep the vertex of the 
parabola right at the q = line, one has to fine-tune at 
least one additional control parameter, besides m. This 
is the case of a tricritical transition p i n ^Si2Si2L2S^. 

In the present work we have restricted ourselves to 
the one-component reactions. It is desirable to extend 
the strategy to reaction-diffusion models with several 
reagents. Each new reagent brings an additional reac- 
tion coordinate and corresponding momentum. E.g., a 
two-component model requires four dimensional phase 
space with d = 3 surfaces of constant energy. The corre- 
sponding classical dynamics may be not integrable, mak- 
ing the phase-space topology rather complicated. The 
situation may be simplified by the presence of conserva- 
tion laws. For example, reaction A + B — > 2B conserves 
the number of particles. This leads to a classical dynam- 
ics with an additional integral of motion, besides energy. 
It is most clearly seen after a canonical transformation 87 
p = e p , q = Qe~ p , which leads to the integral of mo- 
tion Qa + Qb — const. Even with such simplifications 
understanding the full phase-space dynamics of multi- 
component systems remains a challenge. 



arrives at an action of the form 



S = dtd d x 



ip (d t q - DV 2 q 



2 S\ 2 2 

mq + v q J + up q 



(Al) 

We note, in parentheses, that this is the Martin-Siggia- 
Rose action of the following Langevin process with the 
multiplicative noise 77: 



dtq = D\7 2 q + mq 2 — v q 3 + q r}(x, t) 



(A2) 



where {n(x, t)rj(x' , t')) = 2uS(x - x')5(t - t'). 

We shall assume now that the vacuum of the theory 
supports the anomalous average value 



(q 2 (x,t)) = A, 



(A3) 



which is to be determined from the self-consistency con- 
dition. Neglecting the nonlinear fluctuation terms, the 
action acquires the form 

S = Jdt d d x ip (p t - DV 2 + vA) q + uAp 2 - imAp 

(A4) 

With this Gaussian action one can evaluate (q 2 (x,t)) 
and impose the condition 1A3() . This leads to the self- 
consistency equation 



A = uA 



d d k 



rn 



Dk 2 



(A5) 



Without the first term on the RHS, A = (m/v) 2 , this is 
simply the mean-field prediction for (q 2 (x,t)). In d < 2 
this equation has a nontrivial solution even at m = 0. In 
particular at d — 2 one finds 



A (d=2) = 



A 2 D 



, —4ttD/u 



(A6) 
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APPENDIX A: 2CPD MODEL 

Here we consider the 2CPD model (same as PCPD) 
which was shown to exhibit runaway RG flow for d < 2; 
see Sec. IIVBI We shall argue that this behavior of the 
RG indicates a rearrangement of the ground state, such 
that the vacuum supports the anomalous averages of the 
type (g 2 ), similar to the BCS theory. 

In a vicinity of the phase transition the action of the 
model is given by Eq. @ with the reaction Hamiltonian 
(|4UH . To ensure convergence of the functional integral, 
it is convenient to perform rotation of the integration 
contour in the complex p plane: p — > ip. This way one 



where A ~ 1/a is the momentum cutoff. In d < 2 one 
obtains 



A (d<2) 



D 



v \(2-d)D 



2/(2-d) 



(A7) 



As a result, one finds that the "order parameter" A is 
associated with the "coherence length" £, given by 



vA = DC 



(A8) 



Notice that £ is exactly the characteristic spatial scale 
for the breakdown of the RG treatment of Sec. IIVBI see 
Eq. l(4^l) . This consideration suggests that the divergence 
of the RG flow is associated with the development of the 
anomalous average, Eq. 

Let us mention a two-field coupled Langevin descrip- 
tion that has been proposed for the related pair con- 
tact process without diffusion^ and for 2CPD model as 
well2&. In both cases a pair field and a singlet field were 
introduced in order to show a phase transition and esti- 
mate the critical exponents (in fair agreement with those 
reported for microscopic models). More work is needed 
to appreciate the connections (if any) with our anoma- 
lous average microscopic approach. 
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APPENDIX B: FUNCTIONAL 
RENORMALIZATION GROUP 

In some problems one cannot focus on the immedi- 
ate vicinity of the (1,0) point in the phase plane (p,q). 
Instead, one has to keep under consideration the entire 
interval p s [0,1], or even p G [—1,1]. This happens, 
e.g., in parity conserving models, because of the p — > — p, 
q —> —q symmetry. To keep the p interval intact upon 
rcnormalization, one must choose the scaling dimension 
[p] = 0. Since the bare scaling requires [p] + [q] = d, one is 
left with the naive scaling dimension [q] = d. With such 
scaling dimensions one may restrict the reaction Hamil- 
tonian to the lowest powers in q (typically the first and 
second), but one must keep all powers of p. As a re- 
sult, one has to employ the functional RG treatment (cf. 
Ref. [13). 

The generic reaction Hamiltonian for the absorbing- 
state models, mentioned above, is 

H R = f(p)q-g(p)q 2 , (Bl) 

where f(p) = J2n fnP n and g(p) = Y,„ 9nP n are poly- 
nomials of p. From the normalization condition (|13fl it 
follows that 

/(l)=ff(l) = 0. (B2) 

Specific models may possess additional symmetries which 
dictate further restrictions on the polynomials f(p) and 
g(p); e.g., for PC models f(p) is odd, while g{p) is even, 
due to parity conservation: p — > — p, q — > —q. We shall 
keep the presentation general, imposing these additional 
symmetries at a latter stage. Because of the assigned 
bare scaling dimensions [p] = 0, [q] — d, and z = 2, one 
finds [f n ] — 2 and [g n ] = 2 — d. Thus the formal critical 
dimension is d c = 2. 

The one loop renormalization is given by the two- 
vertex loops and leads to the following set of the RG 
equations: 

dlfn = 2/„ - l^Sy] m ( m ~ 1) $n,m+k-2f m gk ,( B 3) 
m,k 

dig n = eg n ~ l^S m(m - 1) S ntm +k-29m9k , (B4) 

where e = 2 — d and S is given by Eq. 145(1 . The fac- 
tors 7}Tn(m — 1) describe the combinatorial number of 
pairs which may form the loop. The S symbols enforce 
the proper number of the external (slow) "legs." Equa- 
tions i|B3|l and i|B4|) may be written as coupled partial dif- 
ferential equations for the functions f(p, I) = J2 n fn{l)p n 
and g(p,l) = Y, n 9n(l)p n - 

dif = 2f-\sgd 2 p f, (B5) 

dig = eg-±Sgd 2 g, (B6) 
For e > 0, Eq. (|B6|) predicts the nontrivial stable fixed- 
point polynomial g*(p), satisfying d 2 g* — 2e/S. In view 



of Eq. I|B2|) the proper solution is 

g*(p) = i(p-l)(p + n), (B7) 

where the parameter k is not specified at this stage. Sub- 
stituting this into Eq. (|B5|I one finds 

d l f = 2f-^(p-l)(p + K)d 2 p f. (B8) 

Since this is a linear equation, one may look for its solu- 
tion in the form 

n 

where the eigenfunctions <j) n (j>) are solutions of the sta- 
tionary equation 

EnMp) = | (P ~ 1)(P + «) d*<t> n (p) . (B10) 

Notice that if one chooses 4> n (j) to be a n-th degree poly- 
nomial, the RHS of Eq. ilBlOfl is also a polynomial of the 
same degree. It is clear then that one may always find a 
solution 4> n (p) as a polynomial of the degree n. This 
means that if one started from a polynomial of some 
degree N, higher powers will not be generated by the 
RG. That is, the sum in Eq. I(B9(I is always confined to 
1 < n < N. 

To find the eigenenergies E n one needs to compare the 
coefficients of the leading power of p on both sides of 
Eq. (|B10l) . This leads to E n = n(n ~ l)e/2. From here 
and Eq. (|B9p one finds that the scaling dimensions of 
the coefficients of the /(p) polynomial at the nontrivial 
fixed-point, Eq. (|B7|I . are 

[f n ]=2-E n = 2-n(n-l)e/2. (Bll) 

The conclusion is that it is sufficient to keep the polyno- 
mials g(p) to be of the second degree [cf. Eq. (|B7(l ] while 
/(p) to be of the lowest possible degree consistent with 
the symmetries of the model. 

In the PC model of Sec. IV Al g(p) must be even and 
thus k = 1, leading to g(p) = 32 (p 2 — 1)- On the other 
hand, /(p) must be odd and at least of degree n = 3 [n = 
1 odd polynomial can not satisfy Eq. (|B2|l ]. Therefore it 
can be chosen to be /(p) = /s(p 3 —p)- At the fixed point 
92 = e /S an d [/s] = 2 — 3e. This justifies treatment of 
Sec. IV Al with the identification v = 32, and u = fy. 

In the reversible model A <-> 2 A of section IVT1 gjjp) = 
v(p 2 — p) and thus k = 0. On the other hand, /(p) = 
u{p 2 — p) which is the proper stationary eigenfunction of 
Eq. IBlOfl : (j>2(p)- As a result, no other terms in f(p) 
polynomial are generated upon renormalization. This is 
consistent with the robustness of the rectangular struc- 
ture. Since u = /2, its scaling dimension according to 
Eq. (EH is [u] = 2 - e. 
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79 In principle, k = 4 may exhibit a marginally non-mean- 
field behavior in d = 1, but we shall not try to investigate 



it here. 

The restricted "fermionic" models may exhibit a different 
critical behavior in low-dimensional systems, than unre- 
stricted models. The most striking example is a restricted 
analogy of the 2CPD model without diffusion (so called pair 
contact process PCP). The PCP model belongs to DP uni- 
versality class (even though there exist an infinitely many 
absorbing statea^S). 

This leaves the possibility that even for d > 2 the RG flows 
towards a strong-coupling fixed point, if the bare nonlin- 
earity is sufficiently strong. As a result one may expect to 
observe a non-mean-field behavior for d > 2. 
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